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Extraordinary carrier multiplication gated 
by a picosecond electric field pulse 
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The study of carrier multiplication has become an essential part of many-body physics and 
materials science as this multiplication directly affects nonlinear transport phenomena, and has 
a key role in designing efficient solar cells and electroluminescent emitters and highly sensitive 
photon detectors. Here we show that a l-MVcm~^ electric field of a terahertz pulse, unlike a 
DC bias, can generate a substantial number of electron-hole pairs, forming excitons that emit 
near-infrared luminescence. The bright luminescence associated with carrier multiplication 
suggests that carriers coherently driven by a strong electric field can efficiently gain enough 
kinetic energy to induce a series of impact ionizations that can increase the number of carriers 
by about three orders of magnitude on the picosecond time scale. 
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mpact ionization originates from the Coulomb interaction 
between carriers and is essentially the reverse of the Auger 
recombination process ^ In this process, a highly energetic con- 
duction-band electron (en) colliding with a valence-band electron 
gives two conduction electrons (ei2 + ^22) ^ hole (/z2i)- Th^ ^l^c- 
tron ei i gains energy in the strong electric field, and this energy gain 
implies that carriers driven without being affected by incoherent 
phonon scattering efficiently acquire enough kinetic energy to trig- 
ger a series of impact ionizations, that is, carrier multiplication^'^. 
This multiplication process has an important role in the efficient 
working of photovoltaic nanomaterials^"^, electroluminescent 
emitters^'^ and highly sensitive photon detectors used in optical 
quantum information applications^^. However, despite the fact that 
carrier- initiated impact ionization in a strong electric field critically 
affects non- equilibrium quantum transport phenomena^ the 
elementary process of carrier scattering relevant to ballistic trans- 
port has yet to be clarified. 

Intense terahertz (THz) pulse sources have the potential to pro- 
vide significant insights into ultrafast carrier multiplication through 
purely optical methods, because a THz pulse has a high electric-field 
amplitude ( > 0.1 MVcm~ ^) and lasts a picosecond^^"^^. It has been 
shown that electrons in the conduction band of GaAs undergo bal- 
listic motion without incoherent longitudinal optical (LO) phonon 
scattering for periods on time scales of hundreds of femtoseconds 
i27iC0i^Q~^=ll5fs, where cOiq is the angular frequency of the LO 
phonon) Carriers accelerated ballistically on ultra-short time 
scales can therefore be expected to gain kinetic energy more efficiently 
and thereby trigger carrier multiplication. However, insufficiently 
strong THz pulses generate few carriers^^, and various additional 
phenomena observed by THz and optical absorption measure- 
ments—such as phonon absorption^ ^ intervalley scattering^^ and 
exciton dissociations^^'^^— obscure the carrier multiplication. 

In this article, we highlight the extraordinarily efficient carrier 
multiplication in the typical semiconductor GaAs when it is sub- 
jected to a THz pulse of a strong electric field. We demonstrate for 
the first time that the electric field of a THz pulse, unlike a DC bias, 
can generate a substantial number of excitons that emit near-infra- 
red luminescence without any help from band-to-band photoexci- 
tation. Our results imply that the carriers coherently driven by a 
strong electric field can efficiently gain enough kinetic energy to 
induce a series of impact ionizations that can increase the number 
of carriers by about three orders of magnitude on the picosecond 
time scale. 

Results 

Carrier multiplication with THz pulse excitation. We developed 
an experimental scheme for studying carrier multiplication in 
nominally undoped GaAs quantum wells (QWs) with intense 
THz pulses (Fig. la. Methods). Carrier multipHcation was induced 
by irradiating the QW sample with a nearly half- cycle THz pulse 
lasting 1 ps and having a maximum peak electric field amplitude of 
1.05MVcm~^ (Methods) The THz pulses were focused onto a 
GaAs QW sample that was mounted on the cold finger of a helium- 
flow cryostat. The electric field polarization was perpendicular to 
the stacking direction and along the (100) direction of the sample. 
As the thickness of the QWs (L = 6 [im) is sufficiently thin compared 
with the centre wavelength of THz pulse (-300 |am; Fig. lb), the 
constructive interference inside the sample modifies the strength 
and shape of the incident THz pulse (Methods). 

The field ionization due to the THz pulse caused electrons from 
residual impurity donors to be accelerated along the A line around 
the Brillouin zone centre (F point) in momentum space as shown 
in Figure Ic (see the Methods section for a description of the field 
ionization)^^. When the electric field reached a certain strength, the 
electrons acquired enough energy to induce several impact ioniza- 
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Figure 1 1 Measurement of THz pulse induced luminescence and impact 
ionization process, (a) Generated THz pulses are focused onto the GaAs 
QWs sample, and the luminescence is detected by a CCD camera after 
it has passed through a spectrometer, (b) The geometry of the sample 
interfaces with air (n/\ = 1), QWs (hqw = 3.5) and a quartz substrate 
(r)s = 2.1). We assume here that the QWs with thickness (/. = 6jim) on 
the quartz substrate has a homogeneous refractive index (hqw = 3.5) 
represented by the average of the refractive indices of the wells (n^ = 3.6) 
and barriers (nb = 3.4). CjCt, x) is incident THz electric field from the air. 
(c) Electron-initiated impact ionization transitions in the schematic GaAs 
band structure for momentum in the A direction. The lattice constant a 
of GaAs is 5.6 A, and ±2K/a corresponds to ±1.1xl0^^m~l The diagram 
shows electrons and hole positions before and after the transition at 
the threshold. 



tions and produce substantial numbers of electron-hole (e-h) pairs 
that go on to form excitons around the T point. The impulsive elec- 
tric field of the pulse allowed the unbound e-h pairs to form exci- 
tons that decayed radiatively because the electric field that would 
otherwise ionize the exciton was absent just after the pairs were 
generated^^'^^'^^. The luminescence from the excitons was detected 
by a Hquid-nitrogen-cooled charge -coupled device (CCD) camera 
(Fig. la). 

Excitation intensity dependence. Figure 2a clearly shows that THz 
pulses with different field strengths induce near-infrared lumines- 
cence centred around 1.55 eV from QWs at lOK, even though the 
central photon energy of the pulse (~4meV) is about 390 times 
lower than the luminescence photon energy^^, and the luminescence 
intensity drastically decreases as the electric field e decreases. As 
luminescence is proportional to the created e-h pair density^^, the 
generated carrier density JV can be estimated by comparing the lumi- 
nescence measurement with optical pulse excitation measurements 
(Methods). As shown in Figure 2b, the carrier density N plotted 
as a function of THz pulse fluence / is extremely nonlinear and 
roughly proportional to (oc e^). The minimum and maximum flu- 
ences correspond to electric fields of 0.47 and 1.05 MVcm"^ and 
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Figure 2 | Luminescence induced by intense THz pulses, (a) Spectra 
of luminescence excited at 10 K by THz pulses with peak electric field 
amplitudes of 0.54 MVcm"^ (green dotted line), 0.70 MVcm"^ (blue 
dashed line) and 1.05 MVcm"^ (red solid line). Inset: normalized 
luminescence spectrum excited with the different electric field strengths, 
(b) THz pulse excitation intensity dependence of peak luminescence 
intensity (magenta closed circles). The eye guide (black solid line) 
corresponds to a fourth-power-law intensity dependence. 



the carrier density increases over this range by about three orders 
of magnitude. Moreover, the shape of the emission spectra does not 
have any field dependence (inset of Fig. 2a), which suggests that 
high- density excitation phenomena do not occur at these excitation 
intensity energies^^. 

Origin of luminescence. We attempted to identify the origin of the 
observed luminescence by comparing the luminescences induced 
by a 1.05MVcm~^ THz pulse and by a 3.18-eV optical pulse. As 
shown in Figure 3a, the spectral shapes are almost the same, and the 
peak photon energies are near the heavy-hole (hh) exciton absorp- 
tion peak. The luminescence induced by the THz pulse can there- 
fore be inferred to be from radiative recombination of hh excitons. 
Figure 3b shows the temperature dependence of peak luminescence 
intensities obtained from the THz and optical pulse excitations. The 
temperature dependences are very similar, suggesting that these 
two kinds of pulses induce similar non- radiative decay processes of 
e-h pairs that contribute in almost the same way to their recombi- 
nation processes. The ratio plotted in the inset of Figure 3b shows 
almost no temperature dependence, and this means that the exci- 
ton generation processes driven by THz pulses and optical pulses 
have similar temperature dependences. In contrast to the similar 
temperature dependences, the extremely nonlinear enhancement of 
luminescence intensity as a function of THz field strength 
nevertheless indicates that the e-h pair creation mechanism due 
to THz-pulse irradiation is completely different from the one due 
to photoexcitation. 

Discussion 

We attribute the generation of substantial numbers of carriers to 
highly efficient carrier multiplication triggered by impact ioniza- 
tions rather than to Zener tunnelling. Note that we can exclude 
Zener tunnelling as a cause because the tunnelling occurring under 
a IMVcm"^ DC electric field lasting Ips can account for only a 
small amount of the e-h pairs (10^^ cm~^)^^, far too small to explain 
the observation (Fig. 2b). In the impact ionization process ^ an elec- 
tron in the conduction band {en), with kinetic energy gained from 
the electric field, creates an electron (^12) plus an e-h pair (^22 + ^21) 
(Fig. Ic): 



-^21- 



(1) 



In accordance with energy and momentum conservation, this proc- 
ess is possible only if the electron has a threshold kinetic energy 
higher than the band gap of 1.52 eV. Here, E^h^^^ i'^^e + ^hh)^ 
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Figure 3 | Comparison of luminescence obtained by THz and optical pulse 
excitations, (a) Spectra of luminescence at 10 K excited by a 1.05 MVcm"^ 
THz pulse (red solid line) and optical pulse with centre photon energy of 
3.18eV (grey open circles) (fluence absorbed by the sample is ~6jiJcm~2). 
(b) Temperature dependence of peak luminescence intensities excited 
by THz (red closed circles) and optical pulses (blue open circles). Inset: 
temperature dependence of intensity ratios as calculated by dividing the 
peak luminescence intensity due to a THz pulse by that due to the optical 
pulse (purple open circles). 

(m^ + m}i}i) = 1.7 eV (ref. 30), where ( = 0.067 mg) and m/^/j ( = 0.5 
mo) are the effective masses of the electron and heavy hole, respec- 
tively. As shown in the sketch of Figure 4a, if an electron released 
from an impurity donor can be accelerated to the threshold kinetic 
energy, it can create an e-h pair (Fig. Ic). Subsequently, both the 
electron that lost energy to create the e-h pair and the created e-h 
pair could, in turn, gain energy from the field. This process can con- 
tinue as long as the electric field of the THz pulse remains sufli- 
ciently strong, and so, comparably few impact ionization events can 
nevertheless produce numerous e-h pairs. 

This intuitive consideration of the multiplication dynamics leads 
to a prediction that can be compared with experimental observa- 
tion. As can be seen in equation (1), an electron -initiated impact 
ionization event doubles the number of electrons. That is, given an 
initial electron density ATq, the electron and hole densities after 

(rii) impact ionization events are, respectively, NqX2^^^^ and 

ATq x(2^^^^ -1). We can neglect the contribution of the gener- 
ated heavy holes because their large effective mass would increase 
the threshold energy Eff^ to 2.9 eV. Thus, the increase in impact 
ionizations (Anj) that is due to changing the electric field from 
^min = 0.47 to Cj^ax =1-05 MVcm ~ ^ can be evaluated from the results 
to be log2(iV(emax)/M^min))~9 (red open circles in Fig. 4b). That is, 
at least nine impact ionizations on average occur within the 1-ps 
duration of the THz pulse. Hence, the impact ionization rate for 
the maximum peak electric field is~9/10~^^s~10^^s~^ 

The experimentally obtained impact ionization rate can be 
quantitatively compared with a simple theoretical prediction, 
whereby the motion of an electron in an electric field follows the 
dispersion of the band structure in momentum space. Here, we take 
the number of impact ionization events (ni(£)) to be the number of 
times the average kinetic energy of the electrons reaches the E^h of 
1.7 eV; the GaAs dispersion relation provides this a wavenum- 
ber k of ±2.77x10^ m~ ^ (ref. 26). The change in electron wavenum- 
ber k(t) is solely determined by the electric field £ (t) acting on the 
electrons, because of the negligible influence of phonon scattering 
processes during periods of hundreds of femtoseconds^^: 



dt 



(2) 



where e is the electron charge and h is Planck's constant. Here, we 
used the electric field of a THz pulse with multiple reflections inside 
the sample (Methods). The change in electron wavenumber k{t) 
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Figure 4 | Comparison of experimentally and theoretically obtained carrier densities, (a) The sketch visualizes the distortion in the Coulomb potential 
of donors, causing the potential to widen and free electrons to be released, and subsequent evolution of unbound e-h gas generated by a series of impact 
ionizations into a pure population of excitons emitting luminescence, (b) Electric field dependences of carrier density obtained in the experiment (red 
open circle) and calculations using equation (2) for three different impact ionization rates /i^^'; one is infinity (grey solid line), and the others are derived 
from equation (3) with C = Co = 870ps"^eV"^ (green one-dot-dashed line) and C=5Co (blue two-dot-dashed line). We assumed A/q^IO^-^ cm"-^. For 
the experimental (red open circles) and calculated data (grey solid line), the carrier density N is plotted together with the corresponding (n|). Panels 
(c) and (d) show the electron wavenumber /c(0 calculated by using equation (2) for e = 1.05 and 0.47MVcm"\ respectively. Dashed lines indicate the 
wavenumbers in the range of ±2.77x10^ m~\ where the electron energy corresponds to the threshold energy E'^^ of 1.7 eV determined from the dispersion 
of the GaAs band structure, (e) Normalized electric field of the temporal profile of the incident THz pulse (grey dotted line) and of THz pulse with multiple 
reflections inside the sample (orange solid line). 



during the THz pulse is calculated under the assumption that the 
electrons lose all their kinetic energy when they reach the average 
wavenumber ±2.77xl0^m~^ that is, an infinite impact ionization 
rate of electrons (Fig. 4c-e). 

Accordingly the e-h densities JV(£) after the (ni(e)) impact ion- 
ization events can be estimated to be iVoX2^"^^ (grey solid line in 
Fig. 4b). The theoretical increase in number (Anj) from ^min-0-47 
to £max=l-05MVcm-i is (ni(£niax))-(^l(^min)) = 12, and this 
estimate would also be subject to an uncertainty of around ±30%, 
owing mainly to the error in measuring the strength of the THz 
electric field. Hence, we can say it is close to the experimentally 
obtained (Anj), and it confirms the validity of the experimentally 
derived Ji of -10^^ s~ ^ This high-impact ionization rate means that 
there are too few phonon interactions to prevent carriers from being 
driven coherently and multiplied efficiently during intervals of hun- 
dreds of femtoseconds. This conclusion is consistent with the tem- 
perature independence of the luminescence intensity ratio (see the 
inset of Fig. 3b). If phonon scattering were a significant process on 
the sub -picosecond time scale, the efficiency of carrier generation 
due to THz pulse irradiation would be lower and the ratio would 
decrease with increasing temperature and phonon population. 
The assumption of an infinite /[^^ is reasonable only if the impact 
ionization time Tj = ij^) ~ ^ is sufficiently short compared with the 
time for the electron energy to reach from zero. The similar 
values of (Anj) between the experiment (~9 times) and calculation 

4 



(-12 times) suggest that Tj is sufficiently shorter than T^. If the 
values were otherwise, the experimentally derived (nj) would be 
suppressed and the carrier multiplication could not occur so as to 
induce a 1000-fold increase in carrier density. 

We also theoretically analysed the observed carrier density with 
a finite impact ionization rate which is modelled according to 
the well-established Keldysh formula^^: 

yf=C{E-Ej. E>E,^. (3) 

The corresponding values for electrons are taken from Fischetti 
and Laux^^: C= 870ps~ ^ eV~^ and E^h-l-^eV. As shown in 
Figure 4b, when we use a value of C = 5Cq in equation (3), the cal- 
culation (blue two-dot-dashed line) reproduces the experimentally 
observed carrier density and is close to the curve calculated under the 
assumption of infinite /[^^ (grey soHd line). Here, /[^^ = 2.5x 10^^ s ~ ^ 
can be derived by using equation (3) with C = 5Cq and E= 1.94 eV 
at the maximum first-conduction-band energy between T and X 
points. The value y^ (-10^^ s~^) obtained experimentally is almost 
ten times smaller than the calculated /[^^ The finite time for elec- 
trons to accelerate to can be the reason for the lower experimen- 
tal value. 

The observed saturation behaviour at stronger electric fields in 
Figure 4b can be qualitatively understood to be from an enhance- 
ment in Coulomb scattering among the generated carriers. Figure 5 
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Figure 5 | Results of model calculation. Calculated carrier density N(t) 
as a function of time for different electric field strengths, which can 
be calculated by using Nox2^"'^ and the time development of impact 
ionization events {n\) during the THz pulse; 0.47 MVcm"^ (blue solid line), 
0.70 MVcm"^ (green solid line) and 1.05 MVcm"^ (red solid line). 



plots the carrier density in the QWs sample as a function of time, 
which can be calculated by using ]Vqx2^"^^ and the time develop- 
ment of impact ionization events (nj) during the THz pulse 
(Fig. 4c-e). Under strong electric fields, the substantial numbers 
of e-h pairs generated at the earlier stage of multiplication may 
undergo a significant amount of e-h scattering. The scattering 
decelerates the electrons and reduces their ability to gain kinetic 
energy. Moreover, e-e scattering could lead to a redistribution of 
energy within the electron system, which may in turn reduce the 
average kinetic energy of a given carrier below the threshold and 
suppress the impact ionization events. Eventually, unbound e-h 
pairs generated with an energy conversion efficiency of ~5xl0~^ 
(Methods) may form excitons that have a relaxation time TpL of 
several nanoseconds (TpL as determined by the radiative and non- 
radiative relaxation times; see Fig. 4a). 

In conclusion, we demonstrated extraordinarily high carrier 
multiplication in the typical semiconductor GaAs when it is sub- 
jected to a THz pulse of a strong electric field; the number of carri- 
ers increased by about three orders of magnitude. The carrier mul- 
tiplication increase that occurs with an increase in the THz electric 
field was in agreement with phenomenological theory based on 
an impact ionization model including electron motion in /c- space 
with a pristine band structure. In the future, a full quantum kinetic 
theory treatment of carrier dynamics in a band structure modified 
by an intense THz electric field should be performed to reveal the 
microscopic origin of the carrier multiplication^^. Moreover, an 
experiment with a time-delayed probe pulse would allow for the 
carrier multiplication dynamics to be directly probed. Our find- 
ings of efficient ultrafast carrier multiplication bode well for future 
applications in ultra-high- speed devices such as high-quantum-effi- 
ciency THz-biased avalanche photodiodes that have femtosecond 
resolution and are sensitive to a single photon, and they may also 
be exploited in the development of efficient electroluminescent and 
photovoltaic nanoscale devices. 

Methods 

Samples. A sample of nominally undoped GaAs/Alg^Gao.yAs QWs was grown 
by molecular beam epitaxy on a (lOO)-GaAs substrate. It consisted of 272 periods 
of GaAs QWs (1 1.9 nm) separated from each other by Alo.3Gao.7As barriers 
(10 nm). We removed the GaAs substrate by wet etching and placed the QWs 
film on a Z-cut quartz substrate (Fig. lb). The dominant donor impurities were 
sulphur from an arsenide source. The static electric field Cjon necessary to ionize 
the sulphur donor is estimated as £ion = -E^b/^^B~4kVcm~ ^ where e is the charge of 



the electron and ^nd are the binding energy and Bohr radius of the impurity 
donor in bulk GaAs, respectively. The given by the literature is ~6 meV, and flg 
{ = £QaQmQ/m'^) is estimated to be -lOnm, where £0= 12.4 and m"^ = 0.067 mo are 
the dielectric constant and effective mass of an electron in bulk GaAs, respectively, 
and ao is the Bohr radius in the hydrogen atom^"^. This ionizing field is much lower 
than the minimum peak field of 0.47MVcm~ ^ to see the luminescence. Conse- 
quently with this set up, the donors are completely ionized during the rise time 
of each THz pulse, and the released electrons trigger the carrier multiplication 
process (Fig. 4a). 

Measurements. High-electric-field THz pulses were generated by optical rectifica- 
tion of femtosecond laser pulses in LiNb03 by using the tilted-pump-pulse-front 
scheme^^"^^'^^. For the pump source, we used an amplified Ti:sapphire laser 
that provided a pulse energy of 4 mj, a full width at half-maximum intensity of 
85 fs, a central wavelength of 780 nm and a repetition rate of 1 kHz. The optical 
pulses from the oscillator with a repetition rate of 80 MHz synchronized with the 
amplified pulses were used for electro-optic sampling of THz pulses at the sample 
position^^. The THz pulse energy measured by a pyroelectric detector (MicroTech 
Instruments, Inc.) was 2.3 \i] at an electric field of 1.05 MVcm~ ^ The generated 
luminescence was analysed with a spectrometer with a resolution of about 1 meV 
and detected by a liquid-nitrogen-cooled CCD camera. To prevent stray near-infra- 
red pulse laser light from entering the CCD camera, we put a black polypropylene 
sheet whose Fresnel coefficient for transmission in the THz frequency region was 
0.87 in front of the cryostat window. To change the THz pump fluence, we put 
black polypropylene sheets in the THz-beam path. 

Electric field of the THz pulse Inside the sample. Figure lb shows the geometry 
of the sample interfaces with air {nj^=l) and the quartz substrate (/is = 2.1). We 
assume here that the QWs has a homogeneous refractive index nqy^ that is the 
average of the refractive indices of the wells («^ = 3.6) and barriers (% = 3.4), that 
is, 3.5. The electric field £{t,x) inside the QWs in which multiple reflections take 
place can be written as follows: 

where £i(t, x) is incident THz electric field from the air, iV is the number of round 
trips, c is the speed of light, t =2/(«q^+ 1) is the Fresnel transmittance between 
the air and sample, = inqyq- 1)/(«q^+ 1) is the Fresnel reflectance at the air and 
sample interface, Y2 = (^qw~ ^s)/(^qw + ^s) is the Fresnel reflectance at the sample 
and substrate interface and To = riQy^ Lie is the phase shift. L is the thickness of the 
QWs. Consequently, the ratio of the maximum values of the positive lobe between 
the 6-|im thick and an infinitely thick QWs increases to 1.2. The maximum peak 
electric field inside the sample is estimated to be 0.55 MVcm~ ^ in the case of a 
1.05 MVcm~ ^ electric field incidence. 

Estimation of carrier density. The carrier density N generated by a 1.05 MVcm~ ^ 
THz pulse can be estimated by the following equation, 

CtHz(Wx) 



iV(£ 



^(«eff/2) 



where Cxhz is the total number of carriers and ( = 3.2 |im) is the length of the 
sample's GaAs wells alone (excluding the AlGaAs barriers). The effective diameter 
of the excited spot area agff (-100 |im) is obtained by dividing the spot diameter 
in the electric field a-^^ (-300 |am) by a factor of 2^/2 , as the generated carrier 
density roughly follows an 8th power law of the incident electric field (Fig. 2b). 
We obtained the same PL intensity with a THz pulse with a peak field value of 
^max with an optical pulse with a photon energy p of 3.18 eV, a spot diameter 
^opt of -100|im and absorbed fluence F of 6|ijcm~^. The carrier numbers were 
about the same because the spot diameters of agff and a^pt were almost the same 
value (-100 |im). Accordingly, iV(eniax) = F / epd^ ~4x 10^^ cm~^, where 
e = 1 .6x10" eV ~ ^ J. (The uncertainties of the carrier densities N are typically 
±50%, owing mainly to the errors in the measurement of the optical excitation 
fluence and spot diameter of the THz pulse.) 

Energy conversion efficiency. The energy conversion efficiency Q from the THz 
pulse energy Wj^^^ absorbed in the sample with respect to the impact ioniza- 
tion process Wa can be estimated from the generated carrier density. Part of the 
absorbed energy of the incident THz pulse goes into the impact ionization process; 
the other parts go into free carrier absorption and the transmitted THz pulse. The 
efficiency Q is defined as 



WA 



^max 

);r(fleff/2) rf„ 
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where £th = 1-7 eV is the threshold energy of impact ionization. As the absorbed 
THz pulse energy Wj^^ is almost 1 the conversion energy efficiency Q is 
estimated to be ~5xl0~'*. 
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